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Abstract 

I reconsider the forward-backward asymmetry for flavoured quarks in electron-positron 
annihilation. I suggest an infrared-safe definition of this observable, such that the asymmetry 
may be computed in perturbative QCD with massless quarks. With this definition, the first 
and second order QCD corrections are computed. 



1 Introduction 



The forward-backward asymmetry of bottom and charm production in electron-positron annihi- 
lation is a key observable for precision measurements. It provides a precise determination of the 
effective weak mixing angle sin 2 0^yy, which in turn leads to constraints on the mass of the 
yet unobserved Higgs boson. In a recent report of the electroweak working group [1] on preci- 
sion measurements on the Z resonance, the forward-backward asymmetry for ^-quarks shows the 
largest discrepancy with the Standard Model predictions of about 2.8 standard deviations. This 
motivates a closer re-examination of the forward-backward asymmetries. 

The forward-backward asymmetries for charm and bottom quarks are measured experimen- 
tally with a precision at the per cent level. To match this experimental precision the inclusion 
of QCD corrections in a theoretical calculation is mandatory. The next-to-leading order (NLO) 
corrections to the forward-backward asymmetry have been calculated in [2-6]. The next-to- 
next-to-leading order (NNLO) corrections have been considered by three groups [7-9] in mass- 
less QCD. The last calculation by Catani and Seymour [9] was performed to clarify the dis- 
agreement between the first two calculations. Furthermore these authors pointed out, that the 
forward-backward asymmetry is not infrared-safe, if the direction is defined by the direction of 
the flavoured quark, as it was done in the first two calculations. Nor is the quantity infrared-safe, 
if the direction is defined by the thrust axis, a definition which has been used in the experimental 
analysis [10]. For a fully massive calculation only partial results for the NNLO corrections are 
available at the moment [11,12]. 

In order to reduce uncertainties related to non-perturbative effects, an infrared-safe definition 
of the forward-backward asymmetry is highly desirable. In this letter I reconsider the forward- 
backward asymmetry for flavoured quarks. The first purpose is to suggest a definition of the 
forward-backward asymmetry, which is infrared- safe. This definition is a direct application of 
the infrared-safe definition of a heavy-quark jet, introduced in [13]. An infrared-safe observable 
can be calculated reliably in perturbation theory. The second purpose of this letter is to com- 
pute for the infrared-safe definition of the forward-backward asymmetry the NLO and NNLO 
corrections. It is well known that the genuine two-parton contributions to the corrections can- 
cel in the asymmetry. Therefore the NNLO corrections could be obtained by a dedicated NLO 
calculation, a fact already noted in ref. [7]. In this paper a different and more direct approach 
is taken. The corrections are obtained with the help of the recently developed general program 
for NNLO corrections to the process e + e~ — > 2 jets [14]. It turns out that these corrections are 
small, making the observable suggested here a candidate for high precision measurements at a 
future linear collider. 

This article is organised as follows: Basic definitions related to the forward-backward asym- 
metry are reviewed in the next section and an infrared-safe definition of this observable is pre- 
sented. In sect. 3 a few technical details are outlined for the calculation of the NNLO corrections. 
The numerical results are presented in sect. 4. Finally, sect. 5 contains the conclusions. 
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2 Definitions 



The forward-backward asymmetry is defined by 

A fb = — , (1) 

OS 

where 05 and Oa are the symmetric and antisymmetric cross sections. These are expressed in 
terms of the forward and backward cross sections Op and 05 as follows: 

o s = o F + o B , o a = o f -o b . (2) 

The definition of the forward and backward cross sections Of and Ob requires a little bit more 
care. Loosely speaking, we think about the forward cross section as that part of the total cross 
section in which the flavoured quark is observed in the forward hemisphere, and correspondingly 
we think about the backward cross section as that part in which the flavoured quark is observed 
in the backward hemisphere. One is therefore tempted to define the direction entering the defini- 
tions of Of and Ob as the direction of the flavoured quark. This definition was used in refs. [7-9]. 
An alternative definition defines the direction by the thrust axis, the two-fold ambiguity inher- 
ent in this definition is resolved by singling out the direction which is closer to the flavoured 
quark. Both definitions are not infrared- safe, a fact which was pointed out in ref. [9]. This fail- 
ure of infrared-safeness is related to the q — > qqq splitting function, which occurs in the NNLO 
corrections for the first time. 

To cure the situation requires a more careful definition of the direction entering the determi- 
nation of Op and Og. It should be taken as the jet axis of a flavoured quark jet. Banfi, Salam 
and Zanderighi [13] recently proposed an infrared-safe definition for a flavoured quark jet. Start- 
ing from the DURHAM algorithm [15], they modified the resolution criteria of the DURHAM 
algorithm 

y OURHAM = 2 ( 1 - g C 2 OS9 ^ min(jE 2 ;jE 2 ) (3) 

towards 

flavour 2(l-cos9;-;) j min(E ; 2 ,E?), softer of i, j is flavourless, 
'j Q 2 \ msut(Ef,Ej), softer of i, j is flavoured. 

The resulting jet algorithm is called the flavour-/:^ algorithm. Each pseudo-particle carries an 
additional additive quantum number, which is the flavour number. If two pseudo-particles are 
merged, the flavour numbers are added. The original DURHAM algorithm is infrared unsafe for 
flavoured jets at order (oc 2 ), which can be seen by considering a large-angle soft gluon which 
in turns splits into a flavoured qq pair. The flavour-^ algorithm ensures that these qq pairs are 
merged together first, by not introducing spurious closeness to other particles in the resolution 
variable. 

As a concrete example let us consider the calculation of the forward-backward asymmetry 
for ^-quarks. We associate the flavour number +1 to a &-quark and the flavour number —1 to a 
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&-quark. All other particles (gluons and d-, u-, s- and c-quarks) are assigned the flavour number 
0. If we combine two ^-quarks, the resulting pseudo-particle has flavour number 2, whereas if 
we combine a Z?-quark with an &-quark, the resulting pseudo-particle has flavour number 0. A 
pseudo-particle with a flavour number equal to zero is called flavourless, otherwise it is called 
flavoured. For the calculation of the forward-backward asymmetry for c-quarks, the roles of the 
b- and c-quark are exchanged. 

This brings us to the following definition of the forward-backward asymmetry for flavoured 
quark jets: We first select two-jet events according to the flavour-^ algorithm and a specified 
resolution parameter y cut . The jet axis defines if a jet lies in the forward or backward hemisphere. 
For the forward cross section we require a jet with flavour number greater than zero in the forward 
hemisphere. The backward cross section is defined similar: Here we require a jet with flavour 
number greater than zero in the backward hemisphere. With these definitions of the forward and 
backward cross sections, the forward-backward asymmetry is then given by eq. (1) and eq. (2). 

This definition is infrared-safe and all quantities can be computed in perturbation theory. In 
particular we are interested in the value of the forward-backward asymmetry on the Z-resonance. 
Here the cross section is dominated by the Z-boson exchange. If we neglect the photon exchange, 
the leading order result for the forward-backward asymmetry on the Z-resonance s = m| is given 
by 

( ) ^ 

Apg = —H e JZq, (5) 

where 

2v f a f ,~ 
Af = 2 2 - (6) 

The quantities Vf and a/ are the vector and axial-vector couplings of a fermion to the Z-boson 
and given by 

v f = h-2Qsm 2 Q w , a f = h. (7) 



3 Outline of the calculation 

In this section I give a brief outline of the calculation of the NLO and the NNLO QCD corrections 
to the forward-backward asymmetry. This calculation is based on a general numerical program 
for NNLO (and NLO) corrections to the process e + e~ — > 2 jets [14]. To this order the following 
amplitudes enter: The amplitudes for e + e~ — > qq up to two-loops, the amplitudes — > qgq 
up to one-loop, as well as the Born amplitudes for e + e~ — > qggq and e + e~ — > qqq'q' . As in 
ref. [9] closed triangle diagrams in the O (oc^) -part are neglected. These diagrams are expected 
to give numerically a very small contribution. The program uses the subtraction method to 
handle infrared divergences [16, 17]. For the calculation of the forward-backward asymmetry 
information on the flavour of the final-state partons had to be added. This is straightforward for 
amplitudes involving only one pair of quarks, e.g. for the amplitudes e + e~ — > qq + n gluons. Here 
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I would like to discuss the contribution from the four-quark final-state. The colour decomposition 
of the tree amplitude for — > qqq'q' is 

-5 //flV ^ (oi 2 5 34 - ^o 14 5 32 ) %f (1 , 4, 3, 2) I , 

2,3,4) = co(l)A(°)(l,2,3,4)+co(3)A( (3,4,l,2). (8) 

In this formula 8 ;J denotes a Kronecker delta in colour space. hfi av equals one, if the flavours of 
the quarks q and q' are identical and zero otherwise. The variable co(q) denotes the electroweak 
coupling of the quark q and is given by 

c (q) = -Qi + v e v^ z (s), <Pz{s)= M2 [. rM - (9) 

s - Mf + iY z M z 

v e and v q are the couplings of the Z-boson to an electron and a quark. ( 1 , 2, 3, 4) is the colour- 
ordered partial amplitude. Explicit expressions for these amplitudes can be found for example 
in [18]. In squaring the amplitude it is convenient to split the resulting expression into a leading- 
colour piece and a sub-leading-colour piece and to distinguish the cases where the two pairs of 
quark are identical or not. We therefore write for the squared amplitude, summed over colours, 
helicities and flavours: 



-(0) 



= flf(l,2,3,4)=flf(l,2,3,4)| /c ^ a +^(l,2,3,4)| Zc ^+^(l,2,3,4)| Wiid . 

(10) 



Note that there is no sub-leading-colour contribution for non-identical quarks. The individual 
contributions are given by 

*(l*^4)| fcfBotid = eY^l\ X T(WA)\ 2 , 

M (l G ,2g,3 e ,4g)| /c . d = e V^^(l%i 0) (l,2,3,4)| 2 + |xi 0) (l,4,3,2)| 2 ), 



^(l e ,2g,3 e ,4g)| sx ,. d = eV^^2Re[ % i 0) (l,2,3,4)* % i 0) (l,4,3,2) 



AN 



(ID 



To render this four-particle final-state contribution integrable, we subtract approximation terms, 
which are re-added to the three- and two-parton final- state terms. In general the double-unre- 
solved contribution is given by 

(o) N n ^° = f(o n+2 d^l 2 -o n+1 oda^-o n oda^ (12) 

The notation O o da is a reminder, that in general the approximation is a sum of terms 

Ooda = Y, da (I 3 ) 
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and the mappings used to relate the n + 2 or n + 1 parton configuration to a n + 1 or n parton 
configuration differs in general for each summand. The individual subtraction terms read as 
follows: The NLO subtraction terms are given by 



da 



(o,i) 

3,qqq'q' 
(0,1 



/c.not id 



da 



3,qqq'q' 



da 



(o,i) 

3,qqq'q' 



lc\A 



id 



N 



(0) 



y [ £ 3°(1,3,4) +£ 3 °(2,4,3)] o * 

~ [£?(1,3,4) + E 3 (2,4,3) + E 3 °(3, 1,2) + £ 3 °(4,2,1) 

+£0(1, 3,2) + E 3 °(2,4, 1) + £ 3 °(3, 1,4) + £ 3 °(4,2,3)] o A 
0. 



(o) 



(14) 



Here, JZ^ denotes the Born amplitude for e + e — > qgq and £ 3 denotes the ggg-antenna function, 
which can be found in [17]. The subtraction terms for double unresolved contributions read 



da 



(0,2) 

(0,0), 2,qqq'q' 



, =(ty- l)^fl2(2,4,3,l)o 
/c,not id I 



(0) 



da 



Ic, id 



doc 



(0,2) 

(0,0),2, 9 ^Y 

[5«(2,4,3, 1) + 2*2(4,2, 1,3) +5°(2,4, 1,3) + #2(4,2,3, 1)] 

(0,2) 

(0,0), 2,qqq'q> 

C F 



(0) 



,id 



^ [Cj(2,4,3, 1) +C 4 °(4,2, 1,3) +C«(2,4, 1,3) +C 4 °(4,2,3, 1)] o 



(0) 



(15) 



A^ denotes the Born amplitude for e + e~ — > gg. i?2(<7, <?) and C®(q, q, q, q) denote the four- 
quark double-unresolved antenna functions. Explicit expressions can be found in [14]. Finally, 
the iterated subtraction terms for double unresolved contributions read 



da 



(0,2) 

(0,1), 2,qqq' q 1 



/c,not id 



N 



N z - 1 



(%-l)^[£ 3 (l,3,4)+£ 3 °(2,4,3)]oi^-±A§(l',2',3') ^ 0) (l",2") 



IN 



da 



(0,2) 

(0,1), 2,qqq>q> 



IN 

hid ~42 



[ J B 3 (1,3,4)+ J E 3 (2,4,3)+ J E 3 °(3,1,2)+ J E 3 (4,2,1)+ J E 3 (1,3,2) 



da 



+£ 3 °(2,4,l)+£ 3 °(3,l,4)+i^ 

(0,2) 

(0,1), 2,qqq>q> 



sc, id 



= 0. 



(16) 



A 3 is the NLO qgq antenna function and can be found in [17]. 
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4 Numerical results 



In this section I give numerical results for the forward-backward asymmetry both for ^-quarks 
and c-quarks. The nominal choice of input parameters is N = 3 colours and Nf = 5 massless 
quarks. I take the electromagnetic coupling to be a(m z ) = 1/127.9 and the strong coupling to 
be a s (mz) =0.118. The numerical values of the Z°-mass and width are mz = 91.187 GeV and 
Tz = 2.490 GeV. For the weak mixing angle I use sin Qw = 0.230. The center of mass energy 
is y/s = mz and the renormalization scale is set equal to /j 2 = s. 

We select two-jet events, where the jets are defined by the flavour-^ algorithm. The recom- 
bination prescription is given by the E-scheme. We consider jets with flavour number greater 
than zero. The direction of a jet is given by the jet axis. 

We expand the symmetric cross section 05, the asymmetric cross section c A , the forward 
cross section Of and the backward cross section o B in powers of a s : 

ok = { 1 + ^ B K+(^) 2 C K ^+0(a 3 s ), wherein e {5,A,F,5}. (17) 

Then the forward-backward asymmetry is given by 

A FB = A^^l + ^5 FS +(^) 2 C Ffi ^+o(a3), (18) 



where the leading order result is given by 



o (0) 

and the coefficients B FB and Cpg can be expressed as follows 



(o) 

A? B = \ (19) 



Bfb = B A -B S , 

C FB = C A -C S -(B A -B S )B S . (20) 
The NNLO coefficients can be decomposed into independent colour structures: 

C K = \(N 2 -l)(c K j c +^C Kjnf +^C K , sc y where K e {S,A,F,B,FB}. (21) 

The leading order result for the bottom and charm quark (including the photon contribution) is 
given by 

A { pl b = 0.11161, A [ pl c = 0.08003. (22) 

The leading order result is independent of the jet resolution parameter y cut . The results for 
the QCD corrections to the forward-backward asymmetry of ^-quarks for various values of the 
jet resolution parameter y cut can be found in table 1 . The corresponding results for the QCD 
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ycut 


BFB,b 


^FB,b 


0.01 


-0.070 ±0.005 


— 0.4 ± 0.8 


0.03 


-0.145 ±0.003 


— 1.7 ±0.5 


0.1 


-0.294 ±0.002 


-4.3 ±0.3 


0.3 


-0.512±0.001 


-10.2±0.1 


0.9 


-0.565 ±0.001 


-13.4±0.1 



Table 1: QCD corrections to the forward-backward asymmetry of b-quarks for various values 
of the jet resolution parameter y cut . The NLO correction is given by B F b, while the NNLO 
correction is given by Cfb- 



ycut 


Bfb,c 


Cfb,c 


0.01 


-0.070 ±0.005 


-0.5 ±0.7 


0.03 


-0.145 ±0.003 


-2.1 ±0.5 


0.1 


-0.294 ±0.002 


-4.8 ±0.2 


0.3 


-0.513 ±0.001 


-12.1 ±0.2 


0.9 


-0.565 ±0.001 


-15.9±0.1 



Table 2: QCD corrections to the forward-backward asymmetry of oquarks for various values 
of the jet resolution parameter y cut . The NLO correction is given by B F b, while the NNLO 
correction is given by Cfb- 



ycut 


i(N 2 -l)C F B,b,lc 


jfi(N 2 -l)C F B,b,nf 


4^r(N 2 - l)C FB ,b,sc 


0.01 


-1.3 ±0.8 


0.85±0.11 


-0.20 ±0.09 


0.03 


-3.3 ±0.5 


1.39 ±0.06 


-0.24 ±0.05 


0.1 


-6.6±0.3 


2.07 ±0.04 


-0.06 ±0.02 


0.3 


-11.4±0.1 


0.97 ±0.02 


0.21 ±0.01 


0.9 


-12.7±0.1 


-0.96 ±0.02 


0.32 ±0.01 



Table 3: The results for the NNLO correction to the forward-backward asymmetry of ^-quarks 
for various values of the jet resolution parameter y cut split into the different colour structures. 



ycut 


i(N 2 -l)C FB ,cJc 


4^(/Y 2 - \ )C F B,c,nf 


^{N 2 -\)Cfb^c 


0.01 


-1.3±0.7 


0.77±0.10 


-0.17 ±0.08 


0.03 


-3.5±0.5 


1.26 ±0.06 


-0.26 ±0.06 


0.1 


-6.7 ±0.2 


1.67 ±0.03 


-0.09 ±0.03 


0.3 


-11.7±0.2 


-0.66 ±0.03 


0.21 ±0.01 


0.9 


-12.6±0.1 


-3.53 ±0.03 


0.32±0.01 



Table 4: The results for the NNLO correction to the forward-backward asymmetry of c-quarks 
for various values of the jet resolution parameter y cut split into the different colour structures. 
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corrections to the forward-backward asymmetry of c-quarks are given in table 2. The errors 
are those of the Monte-Carlo integration used in the numerical program and can be reduced by 
increasing the number of evaluations of the integrand. Tables 3 and 4 give the NNLO results for 
the individual colour structures. 

The QCD corrections are small. By comparing the NLO corrections with the ones reported 
in ref. [9] for the quark axis and thrust axis definition, one observes that the corrections to the 
forward-backward asymmetry defined by the jet axis are for all values of y cut smaller than the 
ones reported in ref. [9] . In addition one deduces from the dependence on y cut that the corrections 
are further reduced by enforcing more stringent two-jet cuts. In addition the definition based on 
the jet axis is by construction infrared- safe, as opposed to the definitions by the quark axis or 
the thrust axis. This makes the forward-backward asymmetry defined by the jet axis a suitable 
precision observable, whose QCD corrections are well under control. Therefore this definition is 
a candidate for high precision studies at a future linear collider. 

As a last remark it is worth pointing out that the numerical program used for this letter can 
also be used with slight modifications to obtain the corresponding QCD corrections to the left- 
right asymmetry and the combined left-right forward-backward asymmetry. 

5 Conclusions 

In this article I reconsidered the forward-backward asymmetry for flavoured quarks in electron- 
positron annihilation with a two-fold purpose: The first point was to provide an infrared-safe 
definition of this quantity. The definition is based on the concept of a flavoured quark jet, where 
the direction is defined by the jet axis. With an infrared-safe definition at hand, the observable 
can be computed perturbatively in massless QCD. In the second part of this article I provided 
results for the NLO and NNLO corrections to this observable. These corrections are small, in 
particular if stringent two-jet cuts are enforced. Therefore this definition of the forward-backward 
asymmetry is an ideal candidate for precision studies. 
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